A pHMM-ANN based discriminative approach to promoter identification in prokaryote genomic contexts

The computational approach for identifying promoters on increasingly large genomic sequences has led to many false positives. The biological significance of promoter identification lies in the ability to locate true promoters with and without prior sequence contextual knowledge. Prior approaches to promoter modelling have involved artificial neural networks (ANNs) or hidden Markov models (HMMs), each producing adequate results on small scale identification tasks, i.e. narrow upstream regions. In this work, we present an architecture to support prokaryote promoter identification on large scale genomic sequences, i.e. not limited to narrow upstream regions. The significant contribution involved the hybrid formed via aggregation of the profile HMM with the ANN, via Viterbi scoring optimizations. The benefit obtained using this architecture includes the modelling ability of the profile HMM with the ability of the ANN to associate elements composing the promoter. We present the high effectiveness of the hybrid approach in comparison to profile HMMs and ANNs when used separately. The contribution of Viterbi optimizations is also highlighted for supporting the hybrid architecture in which gains in sensitivity (+0.3), specificity (+0.65) and precision (+0.54) are achieved over existing approaches.

Genomic analyses and cloning of novel chicken natriuretic peptide genes reveal new insights into natriuretic peptide evolution

The natriuretic peptide (NP) family consists of multiple subtypes in teleosts, including atrial, B-type, ventricular, and C-type NPs (ANP, BNP, VNP, CNP-1–4, respectively), but only ANP, BNP, CNP-3, and CNP-4 have been identified in tetrapods. As part of understanding the molecular evolution of NPs in the tetrapod lineage, we identified NP genes in the chicken genome. Previously, only BNP and CNP-3 have been identified in birds, but we characterized two new chicken NP genes by cDNA cloning, synteny and phylogenetic analyses. One gene is an orthologue of CNP-1, which has only ever been reported in teleostei and bichir. The second gene could not be assigned to a particular NP subtype because of high sequence divergence and was named renal NP (RNP) due to its predominant expression in the kidney. CNP-1 mRNA was only detected in brain, while CNP-3 mRNA was expressed in kidney, heart, and brain. In the developing embryo, BNP and RNP transcripts were most abundant 24 h post-fertilization, while CNP mRNA increased in a stage-dependant manner. Synthetic chicken RNP stimulated an increase in cGMP production above basal level in chicken kidney membrane preparations and caused a potent dose-dependant vasodilation of pre-constricted dorsal aortic rings. From conserved chromosomal synteny, we propose that the CNP-4 and ANP genes have been lost in chicken, and that RNP may have evolved from a VNP-like gene. Furthermore, we have demonstrated for the first time that CNP-1 is retained in the tetrapod lineage.

A genomic approach to understanding the cause and effect of annual ryegrass toxicity

Annual Ryegrass Toxicity (ARGT) is a potentially lethal disease affecting livestock grazing on pastures or consuming fodder that include annual ryegrass (Lolium rigidum) contaminated with corynetoxins. The corynetoxins (CTs), among the most lethal toxins produced in nature, are produced by the bacterium Rathayibacter toxicus that uses a nematode vector to attach to and infect the seedheads of L.rigidum. There is little known of the factors that control toxin production. Several studies have speculated that a bacteriophage specific to R.toxicus may be implicated in CT production. We have developed a PCR-based assay to test for both bacterium and phage in ryegrass material and results indicate that there is a correlation between phage and bacterial presence in all toxic ryegrass samples tested so far. This PCR-based technique may ultimately allow for a rapid, high-throughput screening assay to identify potentially toxic pastures and feed in the field. Currently, ~80% of the 45 Kb genome has been sequenced an investigation to further elucidate its potential role in toxin production.Furthermore, specific alterations in gene expression as a result of exposure to CTs or the closely related tunicamycins (TMs), which are commercially available and considered biologically indistinguishable from CTs, will be evaluated for use as biomarkers of exposure. The effects of both toxins will be analysed in vitro using a rat hepatocyte cell line and screened on a low-density DNA micro array “CT-Chip” that contains <100 selected rat hepatic genes. The results are expected to further define the bioequivalence of CTs and TMs and to identify levels of exposure that are related to specific toxic effects or have no adverse effect on livestock.

Familial and genomic analyses of postural changes in systolic and diastolic blood pressure

The physiological adaptation to the erect posture involves integrated neural and cardiovascular responses that might be determined by genetic factors. We examined the familial- and individual-specific components of variance for postural changes in systolic and diastolic blood pressure in 767 volunteer nuclear adult families from the Victorian Family Heart Study. In 274 adult sibling pairs, we made a genome-wide scan using 400 markers for quantitative trait loci linked with the postural changes in systolic and diastolic pressures. Overall, systolic pressure did not change on standing, but there was considerable variation in this phenotype (SD=8.1 mm Hg). Familial analyses revealed that 25% of the variance of change in systolic pressure was attributable to genetic factors. In contrast, diastolic pressure increased by 6.3 mm Hg (SD=7.0 mm Hg) on standing and there was no evidence of contributory genetic factors. Multipoint quantitative genome linkage mapping suggested evidence (Z=3.2) of linkage of the postural change in systolic pressure to chromosome 12 but found no genome-wide evidence of linkage for the change in diastolic pressure. These findings suggest that genetic factors determine whether systolic pressure decreases or increases when one stands, possibly as the result of unidentified alleles on chromosome 12. The genetics of postural changes in systolic blood pressure might reflect the general buffering function of the baroreflex; thereby, the predisposition to sudden decreases or increases in systolic pressure might cause postural hypotension or vessel wall disruption, respectively.

Extensible motif and aberrant nucleotide extraction in genomic sequences

The project focusses on the discovery of conserved DNA sequences in bacterial genomes and comparative analysis of bacterial genomes to elicit evolutionary trends. The outcomes have produced novel techniques for modelling motifs in DNA and the characterisation of evolutionary processes in medically significant bacterial pathogens.

A genomic approach for improved detection of annual ryegrass toxicity

Annual Ryegrass Toxicity is a severe and constant threat to the Australian agricultural industry. Current diagnostic and detection strategies to predict and monitor ARGT are limited. This thesis utilised genomic-based technologies to develop improved strategies for detection of molecular indicators of toxicity in field and livestock to facilitate pre-clinical detection.

Bacterial genomic G + C composition-eliciting environmental adaptation

Bacterial genomes reflect their adaptation strategies through nucleotide usage trends found in their chromosome composition. Bacteria, unlike eukaryotes contain a wide range of genomic G + C. This wide variability may be viewed as a response to environmental adaptation. Two overarching trends are observed across bacterial genomes, the first, correlates genomic G + C to environmental niches and lifestyle, while the other utilizees intra-genomic G + C incongruence to delineate horizontally transferred material. In this review, we focus on the influence of several properties including biochemical, genetic flows, selection biases, and the biochemical-energetic properties shaping genome composition. Outcomes indicate a trend toward high G + C and larger genomes in free-living organisms, as a result of more complex and varied environments (higher chance for horizontal gene transfer). Conversely, nutrient limiting and nutrient poor environments dictate smaller genomes of low GC in attempts to conserve replication expense. Varied processes including translesion repair mechanisms, phage insertion and cytosine degradation has been shown to introduce higher AT in genomic sequences. We conclude the review with an analysis of current bioinformatics tools seeking to elicit compositional variances and highlight the practical implications when using such techniques.

Mutations in the kissing-loop hairpin of human immunodeficiency virus type 1 reduce viral infectivity as well as genomic RNA packaging and dimerization

A stem-loop termed the kissing-loop hairpin is one of the most highly conserved structures within the leader of human immunodeficiency virus type 1 (HIV-1) and chimpanzee immunodeficiency virus genomic RNA. Because it plays a key role in the in vitro dimerization of short HIV-1 RNA transcripts (M. Laughrea and L. Jette, Biochemistry 35:1589-1598, 1996, and references therein; M. Laughrea and L. Jette, Biochemistry 35:9366-9374, 1996, and references therein) and because dimeric RNAs may be preferably encapsidated into the HIV-1 virus, alterations of the kissing-loop hairpin might affect the in vivo dimerization and encapsidation processes. Accordingly, substitution and deletion mutations were introduced into the kissing-loop hairpin of an infectious HIV-1 molecular clone in order to produce viruses by transfection methods. The infectivity of the resulting viruses was decreased by at least 99%, the amount of genomic RNA packaged per virus was decreased by 50 to 75%, and the proportion of dimeric genomic RNA was reduced from >80 to 40 to 50%, but the dissociation temperature of the genomic RNA was unchanged. There is evidence suggesting that the deletion mutations moderately inhibited CAp24 production but had no significant effect on RNA splicing. These results are consistent with the kissing-loop model of HIV-1 RNA dimerization. In fact, because intracellular viral RNAs are probably more concentrated in transfected cells than in cells infected by one virus and because the dimerization and encapsidation processes are concentration dependent, it is likely that much larger dimerization and encapsidation defects would have been manifested within cells infected by no more than one virus.

Proline residues within spacer peptide p1 are important for human immunodeficiency virus type 1 infectivity, protein processing, and genomic RNA dimer stability

The full-length human immunodeficiency virus type 1 (HIV-1) mRNA encodes two precursor polyproteins, Gag and GagProPol. An infrequent ribosomal frameshifting event allows these proteins to be synthesized from the same mRNA in a predetermined ratio of 20 Gag proteins for each GagProPol. The RNA frameshift signal consists of a slippery sequence and a hairpin stem-loop whose thermodynamic stability has been shown in in vitro translation systems to be critical to frameshifting efficiency. In this study we examined the frameshift region of HIV-1, investigating the effects of altering stem-loop stability in the context of the complete viral genome and assessing the role of the Gag spacer peptide p1 and the GagProPol transframe (TF) protein that are encoded in this region. By creating a series of frameshift region mutants that systematically altered the stability of the frameshift stem-loop and the protein sequences of the p1 spacer peptide and TF protein, we have demonstrated the importance of stem-loop thermodynamic stability in frameshifting efficiency and viral infectivity. Multiple changes to the amino acid sequence of p1 resulted in altered protein processing, reduced genomic RNA dimer stability, and abolished viral infectivity. The role of the two highly conserved proline residues in p1 (position 7 and 13) was also investigated. Replacement of the two proline residues by leucines resulted in mutants with altered protein processing and reduced genomic RNA dimer stability that were also noninfectious. The unique ability of proline to confer conformational constraints on a peptide suggests that the correct folding of p1 may be important for viral function.

Service and cloud computing supporting genomic analysis of the mammalian species

 This research focused on building Software as a Service clouds to support mammalian genomic applications such as personalized medicine. Outcomes of this research included a Software as a Service cloud framework, the Uncinus research cloud and novel genomic analysis software. Results have been published in high ranking peer-reviewed international journals.

Specific recognition of the HIV-1 genomic RNA by the Gag precursor

During assembly of HIV-1 particles in infected cells, the viral Pr55(Gag) protein (or Gag precursor) must select the viral genomic RNA (gRNA) from a variety of cellular and viral spliced RNAs. However, there is no consensus on how Pr55(Gag) achieves this selection. Here, by using RNA binding and footprinting assays, we demonstrate that the primary Pr55(Gag) binding site on the gRNA consists of the internal loop and the lower part of stem-loop 1 (SL1), the upper part of which initiates gRNA dimerization. A double regulation ensures specific binding of Pr55(Gag) to the gRNA despite the fact that SL1 is also present in spliced viral RNAs. The region upstream of SL1, which is present in all HIV-1 RNAs, prevents binding to SL1, but this negative effect is counteracted by sequences downstream of SL4, which are unique to the gRNA.

A survey of approaches and frameworks to carry out genomic data analysis on the cloud

High Performance Computing (HPC) clouds have started to change the way how research in science, in particular medicine and genomics (bioinformatics) is being carried out. Researchers who have taken advantage of this technology can process larger amounts of data and speed up scientific discovery. However, most HPC clouds are provided at an Infrastructure as a Service (IaaS) level, users are presented with a set of virtual servers which need to be put together to form HPC environments via time consuming resource management and software configuration tasks, which make them practically unusable by discipline, non-computing specialists. In response, there is a new trend to expose cloud applications as services to simplify access and execution on clouds. This paper firstly examines commonly used cloud-based genomic analysis services (Tuxedo Suite, Galaxy and Cloud Bio Linux). As a follow up, we propose two new solutions (HPCaaS and Uncinus), which aim to automate aspects of the service development and deployment process. By comparing and contrasting these five solutions, we identify key mechanisms of service creation, execution and access that are required to support genomic research on the SaaS cloud, in particular by discipline specialists. © 2014 IEEE.

Selected approaches and frameworks to carry out genomic data provision and analysis on the cloud

While High Performance Computing clouds allow researchers to process large amounts of genomic data, complex resource and software configuration tasks must be carried out beforehand. The current trend exposes applications and data as services, simplifying access to clouds. This paper examines commonly used cloud-based genomic analysis services, introduces the approach of exposing data as services and proposes two new solutions (HPCaaS and Uncinus) which aim to automate service development, deployment process and data provision. By comparing and contrasting these solutions, we identify key mechanisms of service creation, execution and data access required to support non-computing specialists employing clouds.